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Abstract 

In the double beta decay experiment NEMO 3 a precise knowledge of the background 
in the signal region is of outstanding importance. This article presents the methods 
used in NEMO 3 to evaluate the backgrounds resulting from most if not all possible 
origins. It also illustrates the power of the combined tracking-calorimetry technique 
used in the experiment. 
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1 Introduction 



NEMO 3 is a currently running experiment located in the Laboratoire Souter- 
rain de Modane (LSM) searching for neutrinoless double beta decay {(3(3Qv). 
This decay is a uniquely sensitive probe of the mass and charge conjugation 
properties of the neutrino. Should it be observed, it would demonstrate that 
at least one neutrino is a massive Majorana particle. NEMO 3 is also able to 
detect the rare second order weak double beta decay with its two accompa- 
nying neutrinos (/?/52z^), as well as the non-standard neutrinoless decay with 
Majoron emission {[3j3x)- The three decay modes are distinguishable experi- 
mentally by the energy sum distribution of the two beta particles. Although 
the observation of neutrinoless double beta decay is the goal of NEMO 3, its 
aim is also to measure (or give limits on) half-lives of the other double beta 
decay processes. 

The NEMO 3 detector provides the direct detection of two electrons from the 
decay by the use of a tracking device and a calorimeter (see Fig. 1). It has 
three integrated components: 

- a foil consisting of different sources of double beta emitters {^°'^Mo (6914 g), 
s^Se (932 g), ii^cd (405 g), ^^oXe (454 g), natural Te (614 g of TeOz), ^^°Nd 
(37 g), ^^Zr (9 g), ^^Ca (7 g)} and pure copper (621 g); 

- a tracking volume based on open Geiger cells; 

- a calorimeter made of plastic scintillator blocks with photomultiplier tubes. 

A magnetic field created by a solenoidal coil surrounding the detector pro- 
vides identification of electrons by the curvature of their tracks. Besides the 
electron and photon identification, the calorimeter measures the energy and 
the arrival time of these particles while the tracking chamber can measure the 
time of delayed tracks associated with the initial event for up to 700 /is. A full 
description of the NEMO 3 detector and its performance can be found in [1]. 
The first results were pubhshed in [2] -[4]. 

The most significant concern in double beta decay experiments is the back- 
ground. Although the candidates for double beta decay are selected as events 
with two electron (2e~) tracks with a common origin on the source foil, there 
are certain non (313 processes that can mimic the 2e~ topology. According to 
its origin, the background in NEMO 3 is divided into two categories: 

- the "internal background" having its origin inside the double beta decay 
source foil; 

- the " external background" coming from all radioactive sources located out- 
side of the foil. 

The "internal background" is mainly due to the presence of radioactive iso- 



3 




Fig. 1. Schematical view of the NEMO 3 detector. 

topes from the 238]j g^^^j ^^^Th decay chains. The dominant mechanisms leading 
to the 2e^ topology are (see Fig. 2): 

- f3 decay accompanied by an internal conversion; 

- P decay followed by M0ller scattering; 

- /3-7 cascades in which a 7-ray undergoes Compton scattering. 

The first mechanism is the dominant one. Particularly troublesome are the iso- 
topes with large values such as ^o^Tl (Q;3=4.99 MeV) and ^^^Bi (Q^=3.27 MeV). 
Fortunately, for these isotopes there are good estimates of the activities ob- 
tained from identifiable topologies {ea for ^^^Bi, 677 and 6777 for ^°^T1). 
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Great care was taken in the production and subsequent purification of the en- 
riched materials, as well as diiring the source foil production and the mounting 
of foils in the detector so as to keep any contamination to a minimum given 
the strict radioactivity limits. 

In the region where a signal of neutrinoless Pl3 decay is expected, the allowed 
I3(32p decay can be an important fraction of the background. Its contribution 
depends upon the j3l32v half-life and the energy resolution of the detector. It is 
therefore important to carefully study the background in a large energy region 
where the /?/?2z/ decay takes place in order to obtain a good measurement of 
the half-hfe of this process. 

A component of the external background producing events similar to the in- 
ternal background is caused by the presence of radon and thoron inside the 
detector. These elements are highly diffusive radioactive gases. They are out- 
gased in the air from the rock walls of the experimental hall and can enter 
the detector either through tiny gaps between sectors or through gas pipe 
joints. The progeny of radon and thoron produces 7-rays and /3-decays ac- 
companied by internal conversion (IC), M0ller or Compton scattering. If such 
an event occurs on or near a foil and appears with a 2e~ topology it becomes 
indistinguishable from a double beta decay candidate. 






beta + IC 



beta + M0ller 



beta + Compton 



• = radioisotope (3 = electron from beta decay IC = interna! conversion 

Fig. 2. Internal background production in the source foil. 






pair creation 



Compton + Compton 
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Compton + M 0ller 



Fig. 3. External background production in the source foil. 



Another component of this background is due to external 7-rays interacting 
inside the foil. These 7-rays are of different origins: 
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Fig. 4. Decay chain of the radioactive family of ^^^U. The half-hves and decay 
energies are taken from [5]. 

- 7-rays inside the laboratory, mostly coming from the rock walls; 

- neutron interactions in the shield and material of the detector; 

- radioactive isotopes present in the detector materials despite the rigorous 
selection done during the detector construction; 

- presence of radon in the air surrounding the detector. 

The interaction of 7-rays in the foil can appear like 2e~ events by e^e~ pair 
creation with misidentification of the charge, double Compton scattering or 
Compton scattering followed by M0ller scattering (see Fig. 3). 

In this article the methods used to evaluate the various backgrounds are pre- 
sented. The very pure copper foils (OFHC) are used to prove their validity. 
The experimental data of the NEMO 3 detector were used to perform the 
background measurements. 



2 Natural radioactivity inside the tracking volume 



2.1 ^'^^Rn measurement inside the tracking chamber 



The most bothersome external background comes from radon. As shown in 
Fig. 4, 2i4Bi is one of the descendents of '^^'^Rn. The /^'decay of ^i^Bi to ^i^Po 
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is generally accompanied by several photons which can mimic a (3(30u event 
given its large (3,g value. The ^^"^Po has a half-life of 164 //s and it disintegrates 
to ^^°Pb via a-decay. 

The ejected alpha particle from the decay of ^^^Rn can free several electrons 
from the ^^^Po atom transforming it to a positively charged ion. Diffusing 
through the gas of the tracking chamber this ion may be neutralized by dif- 
ferent processes, such as recombining with negative ions in the gas, a charge 
transfer by neutral molecules with a small ionization potential or the cap- 
ture of an electron created during the gas discharge near the open Geiger 
cell wires [6,7]. It is difficult to predict the proportions of neutral and charged 
atoms of ^^^Po in the gas of the NEMO 3 tracking chamber. However, following 
some earlier studies [8,9,10] one can suppose that the proportion of neutral 
atoms of ^^^Po which are in the gas is small and that the majority of the 
charged ^^®Po is deposited on the surfaces of the cathode wires of the Geiger 
cells. Here the measurement of the radon activity in the tracking chamber is 
done under the assumption that the descendents of ^^^Rn are deposited on the 
wires. It should be noted that the final background resulting from the presence 
of ^^^Rn does not depend critically upon this assumption. If one supposes that 
the descendants of this gas are uniformly distributed in the tracking chamber 
the background estimation remains unchanched within the errors. A possi- 
ble deposition on the foil surfaces is also taken into account in the analysis 
performed to measure internal foil contamination by ^^^Bi (see section 4). 

2.1.1 Event selection 

In the NEMO 3 experiment, the data acquisition allows one to readout signals 
from delayed tracks in order to tag the a particles from ^^^Po decays associated 
with electrons from ^^^Bi decays. The ^^^Bi decays followed by ^^^Po decays 
(hereafter referred to as BiPo events) are used for the measurement of the 
radon activity. Typical examples of such events with detected electron and 
cu-particle tracks are shown in Fig. 5. Two types of spurious events can appear 
in a sample of BiPo candidates. They are due to: 

- random coincidence of two independent events closely localized in space and 
occurring inside the 700 fxs time window. The delay time distribution is flat 
in this case. 

- a single event accompanied by one or more delayed signals caused by refiring 
of neighbouring Geiger cells. The number of this kind of event decreases 
sharply with the delay time, so their distribution is nearly exponential. 

In order to reduce the contribution from refirings, the following cuts are ap- 
plied: 

- for events with only one delayed signal, the delay must be greater than 
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Fig. 5. Examples of BiPo event candidates, viewed from the top. In each event 
the prompt track is shown in red (electron), delayed track - in green (a-particle). 
The two tracks have a common vertex: a) on the source foil, b) inside the tracking 
volume. Scintillator blocks located at top and bottom (wall) are shown as magenta 
(blue) boxes. 

90 fis; 

- for events with more than one delayed signal, the delay must be greater 
than 30 fis. 

It is required that the delayed signals have to be located close to the "prompt" 
(i.e. electron) track or to the event vertex on the source foil. The group of 
selected delayed signals must be within 2.1 //s (corresponding to the maximum 
transversal drift time in the Geiger cells) and follow a straight line in space in 
order to exclude random coincidences of delayed events of different origins. In 
the tracking chamber the maximum range of the 7.7 MeV a-particle emitted 
by ^^■^Po is 36 cm. This value is also taken into account in the selection of 
events. 

Applying these criteria the mean efficiency to select a BiPo event produced 
on a wire surface has been estimated by a Monte Carlo (MC) simulation to 
be 16.5%. In this work all simulations [12] are based on GEANT3 [11] using 
DECAYO [13] as event generator. 

The time distribution of the delayed tracks (see Fig. 6) provides an efficient 
way to validate the quality of the event selection. The fit to the distributions 
allows the half-life of ^^^Po to be evaluated. The proportions of events due to 
refirings and random coincidences are also determined from the fit. They are 
found to be negligibly small with 0.6±1.3% for the single delayed signal and 
1.9±0.8% for the multiple delayed signals as shown in Fig. 6. In spite of the 
systematic uncertainty of ~1 /is on the result, the half-life of ^^^Po for the 
single delayed signal events is Ti/2(^^''Po) = 162. 9±0.2(stat. only) fis and for 
the multiple delayed signal events is ri/2(^^^Po) = 161. 9±0.8(stat. only) /is, a 
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Fig. 6. The time distributions of events selected with (a) single and 
(b) multiple delayed signals. Each distribution was fit to the function: 
f{time) = PI ■ e"-P3/"i"2 + P2 + P4 • e'^^f^ where PI and P4 are scaling con- 
stants, P2 is the amplitude of random coincidences, P3 the ^^^Po half-life in ^us and 
P5 the time constant of the refirings. 

comparison with the table value of Ti/2(^^^Po) = 164.3±2.0 fis [14] confirms 
that the delayed tracks are due to the a-particles of ^^^Po. 



2.1.2 Rn monitoring 

Using the method described in section 2.1.1 the mean ^^^Rn level in the track- 
ing volume was analysed. The results of the radon monitoring are shown in 
Fig. 7. 
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Fig. 7. The ^^^Rn activity in mBq/m'^ inside the tracking chamber measured on an 
hourly basis. 
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In order to decrease the level of ^^^Rn, a radon reduction factory has been in- 
stalled in the underground laboratory to inject nearly radon free air into a tent 
built around the detector. It became operational at the beginning of October 
2004. Therefore the NEMO 3 data has been divided into two parts according 
to the collection dates. The first part of the data corresponds to acquisition 
done from the beginning of the experiment in February 2003 up to the end 
of September 2004 ("Phase 1"). The second part of the data presented here 
includes the runs from October 2004 up to the end of 2006 ("Phase 2"). Data 
collection continues under the conditions of Phase 2. The mean ^^^Rn level 
was measured to be 37.7±0.1 mBq/m^ for Phase 1, and 6.46±0.02 mBq/m^ 
for Phase 2 inside the tracking volume (only statistical errors are given). 

According to the goal of the analysis the ^^^Rn monitoring provides the pos- 
sibility to select data with different radon levels. 



2.1.3 Spatial distribution of the Bi 

In the preceding paragraph the results of the mean ^^^Bi activity measure- 
ments in the tracking volume were given. For double beta decay studies it is 
important to know the ^^'^Bi distribution close to the source foils. 

The starting point of the electron track closest to the delayed track is used 
to localize the decay point of a BiPo event. The position of such a point is 
defined by the Geiger cell layer number, running from 08 to 00 for the inner 
Geiger cell layers and from 10 to 18 for the outer Geiger cell layers (Fig. 8). 
The vertical position is the measured Z coordinate of the electron's origin. 
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Fig. 8. Lay out of the end-cap showing the Geiger cell layer configuration, top view. 
The inner layers are numerated 00 to 08 starting from the source foil shown in red. 
The outer layers are 10 to 18 as indicated. 



Its azimuthal position is given by the sector number from to 19. 



Origin Z = is in the middle plane of the detector 
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Fig. 9. Plots of the Bi activity as a function of Geiger cell layer, sector number 
and vertical position inside the tracking chamber for data of (a) Phase 1 and (b) 
Phase 2. 
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The 2^^Bi distribut ion inside the tracking chamber for Phase 1 and Phase 2 is 
shown in Fig. 9. For both phases there is a consistent pattern that is scaled 
to the radon activity. Greater activity is observed in layers 03 to 06 and 13 to 
16 and is explained by the large gap between Geiger layers 03-04, 13-14, 05-06 
and 15-16 (Fig. 8). 



In order to detect possible sources of ^^^Rn outgasing or ^^^Bi pollution, the 
spatial distribution of ^^"^Bi in Phase 2 has been compared with Phase 1. The 
results are shown in Fig. 10 and indicate that the distributions of ^^^Bi for the 



A^'^Bi) PHASE 2 / A(^^''Bi) PHASE 1 




12 13 14 15 16 17 18 
Geiger cell layer 



Fig. 10. Ratio of Bi activity for two data samples (a) as a function of Geiger cell 
layer and sector number and (b) vertical position in the Geiger cell layers. 

two phases are similar in the whole tracking volume except at the borders. 
There is a residual ^^^Bi activity near the scintillator walls and end caps. 
The enhanced ^^^Bi activity near the foil extremities (especially in sector 3 
with molybdenum foils) probably originates from the foil holders and/or the 
scintillator surfaces. In order to exclude the top and bottom parts with the 
"extra" ^-"^^Bi activity, only events with vertices of |Z| < 120 cm are taken into 
account in the following. 

In the Pf3 analysis most of the background due to ^^^Rn appears to come from 
the regions close to the source foils. Therefore the observed non-uniformity of 
^^^Bi distribution along the vertical direction for distant Geiger cell layers is 
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Table 1 

The ^^^Bi activity in mBq for the Geiger cell layer closest to the copper foils of 
sector 0. The average activity for the 20 sectors is also presented. Statistical errors 
are given. 





A(2i4Bi),mBq 


Phase 1 


Phase 2 


Inner foil side, sector 
Outer foil side, sector 


724±10 
598±8 


134±4 

loizta 


20 sectors average 


700±1 


140±1 



neglected. While the variations in the azimuthal direction from scctor-to-scctor 
are accounted for. The mean activities measured for each sector layer-by-layer 
are used. The test of the background model (see section 5) is performed with 
copper foils located in sector 0. For this sector the results of the ^^"^Bi activity 
measurements for the layers closest to the foil are given in Table 1. 

2.1.4 ^^^Rn activity measurement using 67 events 

It is also possible to detect ^^''Bi using events. A large fraction of the ^^''Bi 
decays is accompanied by a high energy 7-ray. However, the 67 events are 
contaminated by external 7-rays which Compton scatter on the wires of the 
Geiger cells (see section 3). 

To suppress this background contribution, only events with a 7-ray of energy 
greater than 1 MeV are used. In order to select events originating from the 
tracking volume and not from the source foils, only electrons with their starting 
point on the Geiger cell layers 01, 11, 02 and 12 are analyzed. The results for 
the mean ^^^Bi activity for the second and third Geiger cell layers using ea 
and ej topologies are shown in Table 2. 

Table 2 

Results of the ^^^Rn measurements on the second and third Geiger cell layers of the 
tracking chamber using two event topologies are shown. Statistical errors are given. 







A(214 


Bi), mBq 




Geiger cell layer 


02 


01 


11 


12 


ea 


598±3 


701±3 


706±3 


800±3 


67 


688±7 


624±7 


645±5 


727±6 



The method involving 67 events has a larger background with approximately 
three times smaller detection efficiency compared to the method using the 
delayed tracks. However, the 67 events allows one to address the measurements 
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of delayed tracks and to estimate the systematic error on ^^^Rn activity inside 
the tracking chamber to within 10%. 

2.2 '^'^^Rn (^^^Tl) activity measurements inside the tracking chamber 

If ^^°Rn is present in the gas of the tracking chamber, it constitutes a source of 
^°^TL Given the high Qp value, "^^^Tl is a serious concern for neutrinoless dou- 
ble beta decay search. The beta decay of ^ns^l is almost always accompanied 
by a 7-ray of 2615 keV from the first excited state of ^"^^Pb [5]. In the case of 
the de-excitation by IC electron emission one can observe two electrons with 
a total energy of approximately 3 MeV which can mimic ■'^'^^Mo and ^^Se neu- 
trinoless double beta decay. Therefore the thoron content of the gas mixture 
has been studied by means of the ^os^l beta decay detection. 

2.2.1 Selection criteria 

The beta decay of ^os^l is mainly accompanied by two or three 7-rays. There- 
fore 677 ^-iid 6777 topologies are used in the analysis. The event vertex is 
defined by the origin of the electron. In order to reject events coming from 
the source foil only events with vertices on the Geiger cell layers 01—04 and 
11—14 are analyzed. Tracks starting in planes 05 to 08 and 15 to 18 are too 
short to allow accurate event selection by time-of- flight. 

For both topologies the electron has an energy greater than 200 keV and each 
7-ray an energy greater than 150 keV. The time-of-flight method is used to 
ensure the common origin of all the particles involved. In order to reduce 
backgrounds the 7-ray with the highest measured energy (E^) is required to 
have > 1700 keV and the energy of the electron must satisfy the condition 
Ee > (4200 - EE^) keV. 

2.2.2 Results of the measurement 

The data has been analyzed separately for the two phases and the two topolo- 
gies. The results of the measurements are shown in Table 3. There is good 
agreement between the results obtained for the 677 and 6777 topologies. 
The average activities resulting from them are 3.5±0.4 mBq for Phase 1 and 
2.9±0.4 mBq for Phase 2. The systematic uncertainty is estimated to be less 
than 10% and comes mostly from the 7-ray detection efficiency. A compar- 
ison of the results for the two phases shows no strong difference due to the 
detector environment but indicates a possible weak outgasing of thoron inside 
the tracking volume. Taking into account the 35.94% branching fraction of 
^°^T1 in the chain, the measured '^^^Tl activity of ~3 mBq corresponds 
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Table 3 

The number of observed events, number of estimated background events, signal 
efficiency and tfie results of the measurements of ^^^Tl activity inside the tracking 
chamber found for both phases via 677 and 6777 topologies. 



Topology 


N (observed) 


N (estimated bgr) 


Eflf, % 


A(208t1), mBq 


Phase 1 


677 


342 


22.4 


0.26 


3.5±0.4 


6777 


63 


1.8 


0.05 


3.3±0.5 


Phase 2 


677 


322 


6.6 


0.24 


2.8±0.3 


6777 


79 


1.2 


0.05 


3.5±0.5 



to ~8 mBq of ^^'^Rn. MC simulations show that this low level of thoron inside 
the tracking device is much less of a concern for neutrinoless double beta decay 
than the presence of radon. For both phases the background from thoron in the 
2e~ channel is more than one order of magnitude lower than the background 
originating from radon. 

2.3 "^^^Ph (^^^Bi) activity inside the tracking chamber 

One probable source of background with its origin on the Geiger cell wires is 
the /?-decay of ^^^Bi from ^lopb (T1/2 = 22.3 y) in the ^^su decay chain (see 
Fig. 4). Given its Q/3 of 1.16 McV, this radioactive isotope is of no concern for 
this neutrinoless double beta decay search, but it must be considered in the 
precise measurement of the two neutrino double beta decay spectra. 

One electron events with an energy greater than 600 keV and their vertices 
associated with Geiger cells arc selected to measure the ^^"^Bi activity on the 
wire surfaces. The results for both phases are in Fig. 11 and Table 4. 



Phase 1 Phase 2 a (^'°Bi), Bq 




Fig. 11. Plot of Bi activity in Bq on the Geiger cell wire surfaces as a function 
of the event vertex position for (a) Phase 1 and (b) Phase 2. 
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Table 4 

Activity of ^^'^Bi in Bq on wire surfaces of Geiger cell layers which are next to the 
foils. 



Sector 





1 


2 


3 


4 


Phase 1 


6.21±0.37 


2.20±0.11 


2.95±0.28 


3.52±0.21 


2.60±0.17 


Phase 2 


5.83±0.35 


2.10±0.12 


2.76±0.25 


3.31±0.20 


2.41±0.13 


Sector 


5 


6 


7 


8 


9 


Phase 1 


12.62±0.73 


3.26±0.17 


5.09±0.52 


3.29±0.25 


1.89±0.24 


Phase2 


11.84±0.63 


3.17±0.16 


4.71±0.42 


3.13±0.19 


1.85±0.14 


Sector 


10 


11 


12 


13 


14 


Phase 1 


2.86±0.23 


2.10±0.24 


4.51±0.49 


16.72±0.93 


2.21±0.14 


Phase 2 


2.71±0.18 


2.00±0.13 


4.28±0.37 


17.00±1.16 


2.11±0.14 


Sector 


15 


16 


17 


18 


19 


Phase 1 


12.91±0.70 


18.37±0.99 


14.48±0.73 


2.48±0.12 


3.87±0.44 


Phase 2 


12.17±0.63 


17.43±0.93 


13.80±0.70 


2.40±0.12 


3.84±0.41 



Approximately the same activity values and spatial distribution are measured 
in both data samples, while a large variation of ^■^'^Bi activity from one sector 
to another is observed. The origin of the non-uniformity in ^^°Pb deposition 
on the wires is most probably due to the different histories of the wires and 
conditions during the wiring of the sectors. 



3 External 7-ray flux 

The external 7-ray flux is one of the sources of 2e~ events and therefore a 
background for double beta decay. With the NEMO 3 data it is possible to 
measure this flux using the events resulting from single or double Compton 
scattering. 

If an incoming 7-ray undergoes Compton scattering inside a scintillator block 
leaving the scattered electron unseen in the tracking chamber and subsequently 
rescatters in the foil ejecting an electron which hits a scintillator block, an 
event of 67 topology is observed, see Fig. 12a. Such an event from a double 
Compton scattering is identified as an " e7-external" event as opposed to an 
" 67- internal" event with both particles coming from a decay or an interaction 
produced inside the foil. Using time-of-flight measurements and the timing 
properties of the detector it is easy to distinguish an e7-external event from 
an e7-internal event. 

When a single Compton scattering of the incoming 7-ray occurs very close 
to the surface of the scintillator block, the scattered electron, provided it has 
sufficient energy, crosses the tracking chamber. The topology of such an event 
is a pair of tracks of opposite curvature sign coming from the foil. Using 
the track curvatures and time-of-fiight information it is easy to select such 
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Fig. 12. Examples of (a) e7-external and (b) crossing electron events. Presumed 
7-rays shown by wavy lines are superimposed to the event display. 

"crossing electron" events, see Fig. 12b. A possible background for this type 
of event is from /3-emitters at the detector's inner surface. 



3. 1 'f-ray flux from surrounding rocks 

In the LSM underground laboratory there is a significant flux of 7-rays com- 
ing from natural radioactivity in the surrounding rocks. This mostly includes 
gamma radiation from ^°K, ^^"^Bi and ^os^l decays. A passive shield, surround- 
ing the detector, has been constructed to reduce this source of background. 
The shield consists of low radioactivity iron plates and tanks filled with bo- 
rated water, see Fig. 1. 

The 7-ray energy spectrum in the LSM laboratory has been measured using 
a NaT detector. In Table 3 of reference [15] one can find intensities of the 
2.61 MeV and 1.46 MeV 7-ray lines of ^os^l and decays, respectively. The 
intensities of the lines corresponding to ^^^Bi decay have not been evaluated 
in that work, however, with the published spectrum, one can estimate the 
strength of the 1.76 MeV 7-ray line and use it as a reference. Relative ratios 
between different ^^^Bi lines can be found in an earlier work [16]. Although 
these measurements were done for the Gran Sasso underground laboratory, 
one can assume that the rock composition of the two sites is similar, and 
therefore the attenuation is similar too. The result of this compilation for the 
7-ray flux is summarized in Table 5. 

In order to evaluate the corresponding background in NEMO 3 a two stage 
MC simulation has been done. First, a simple MC simulation was performed 
to estimate the attenuation of the NEMO 3 shielding (iron plates and water 
tanks) for the spectrum from Table 5. The MC provided energy and angular 
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Table 5 

Simplified model of 7-ray flux in the LSM underground laboratory. 



Isotope 


T-rav enerffv. keV 


Flux, cm ^ 


40k 


1461 


0.1 


208 J 


2615 


0.04 


2143- 


1764 


0.05 




1600 


0.026 




1300 


0.041 




1120 


0.046 




609 


0.109 


Total 




0.411 



spectra ( dcos(6>) ) 7-rays after the shielding with a total attenuation 

factor of 3.5 ■ 10-^ 

The energy and angular spectra were then parametrized for a MC simula- 
tion with all the components of NEMO 3. This simulation showed that, after 
cuts were applied, the 7-ray flux from the laboratory accounts for ~2% of 
the total measured external background, see Fig. 13. Consequently it can be 
neglected and taken into account by a slight adjustment of other components 
in the external background model. The most important contribution expected 
due to the PMT radioactivity measured with the HPGe detectors [1] is also 
demonstrated in this figure. 




ETOT(MeV) E^(MeV) 

Fig. 13. The sum of the 7-ray and electron energy (Etot), and the electron energy 
(Eg) of the e7-external events in Phase 2 data. The expected contributions due only 
to the 7-ray flux from the laboratory (see Table 5) and from the photomultiplier 
tubes are also shown. 
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3.2 Neutrons 



Neutrons can contribute to the external background via the neutron capture 
process resuhing in emission of 7-rays. The neutron flux in the LSM has been 
measured [17], [18] and originates from spontaneous fission and (a,n) reactions 
due to trace amounts of uranium in the rocks. The NEMO 3 neutron shield 
thermalizes fast neutrons with energies of a few MeV and suppresses thermal 
and epithermal neutrons [1]. 

A series of calibration runs with an AmBe neutron source in the vicinity of 
the detector has been done to check the shield's efficiency. These runs may 
also be used to evaluate the neutron background for /3/?2z/ measurements. The 
energy sum distribution of crossing electron events for Phase 1 data is shown 
in Fig. 14. The distribution for a run with the AmBe source is superimposed 
on it. A pronounced tail at energies up to 8 MeV is seen in both distributions 
and is a characteristic feature of neutron captures in the detector walls. This 
tail above 4.5 MeV is used for normalization of the neutron source distribution 
(Fig. 14). Attributing all the data in this energy region to neutrons one can 
see that the neutron contribution to the low energy part of the spectrum 
measured without the AmBe source is very small. In this way one finds that 
the neutron background can amount to 0.03% of the total at energies below 
4 MeV. Therefore this background can be neglected for p(32i> analysis. 




0123456789 10 



ETOT(MeV) 

Fig. 14. Energy sum distribution of crossing electron events for 404 days of data 
(Phase 1). The superimposed histogram represents the distribution obtained for 25 
hours of data with an AmBe neutron source and scaled by the factor 6.9 • 10^^. 

3.3 External j-ray flux model 

The principal source of the external background is the natural radioactivity 
of the detector components. The dominant one is due to the PMTs contami- 
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nation by ^^^Ra, ^^^Ra and as known from the results of HPGe detector 
measurements [1]. It is addressed in the MC simulation with decays of ^^^Bi, 
208T1, 228Ac and inside the PMT glass. However, the use of the PMT 
activities allows one to reproduce roughly only half of the observed experi- 
mental e7-external events, see Fig. 13. The small fraction of each component 
of NEMO 3 was measured extrapolating the results to the total amount of 
the material. Therefore these activities are used as free parameters in order to 
fit the experimental data. The presence of these isotopes and of ^°Co in other 
parts of the detector have also been taken into account. The hst of potential 
isotopes providing the low energy 7-rays is long. It includes for example ^^Mn, 
^^Co, ^^Zn, ^^^Cs, and ^^^'"Pa. In order to take them into account and to im- 
prove the description of the low energy part of the gamma spectrum, a flux of 
1 MeV 7-rays was simulated at the external surface of the calorimeter. 

Not all components of the apparatus or sources of background have been 
considered so the results of the fit described below should not be interpreted 
as a measurement of the internal contamination of the corresponding elements 
of the detector. The purpose of this study is to provide an effective model able 
to reproduce the external 7-ray flux. 

Table 6 

Components of the external background model. 



Components 
of NEMO 3 


Activity (Bq) 


40k 


214Bi 


208 fpj 


6OC0 


Photomultiplier tubes 


1078. ± 32. 


324. ± 1. 


27.0 ± 0.6 




Plastic scintillators 


21.5 ±0.9 








/x-metal PMT shield 








14.6 ±2.6 


Iron petals 


100. ± 4. 


9.1 ± 1.0 


3.1 ±0.5 


6.1 ± 1.8 


Copper on petals 








47.6 ± 7.8 


Internal tower 








18.4 ± 0.8 


Iron shield 




7360. ± 200. 


484. ± 24. 





The external background was evaluated using two types of events, e7-external 
and crossing electrons. A global fit of the following distributions was per- 
formed: the total energy released in the calorimeter, the energy deposited by 
the electron, the cosine of the angle between the directions of the electron and 
7-ray for the e7-external events, the total energy and finally the energy of the 
electron after crossing the tracking volume for the crossing electron events. In 
order to describe the external 7-ray fiux coming into the detector in more de- 
tail three sets of histograms were produced. These are for the incoming 7-ray 
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signals detected by a counter at the internal wall, external wall and end-caps. 

First the external background model parameters are fixed from the Phase 2 
data where one can neglect the radon in the air surrounding the detector 
given its very low level (~ 0.1 Bq/m^). The activities of the background model 
presented in Table 6 and the 1 MeV 7-ray fiux of 0.446 m~^s~^ refiect the data 
for the incoming 7-rays detected by the calorimeter of NEMO 3. In Fig. 15 
the energy sum (ii^TOT) distributions of the e7-external events are shown. 
The crossing electron channel requires an additional source of electrons that 
was found to be possible to reproduce with 7.3 mBq/m^ of ^^^'"Pa and 120 
mBq/m^ of ^^°Bi on the surface of the scintillator blocks. The results of the 
fit for the whole detector for Phase 2 data are presented in Fig. 16. In the 
Fig. 16(a) one can see the distribution of the total energy measured in the 
calorimeter for crossing electron events, where the beta emitters provide a 
noticeable contribution. However, they are negligible in the case of e7-external 
events for which one can see the distributions of the energy sum of the electron 
and 7-ray in Fig. 16(b) and of the detected electron energy in Fig. 16(c). 

Then the Phase 1 data were fit with the model including the contribution 
of radon and thoron in the air thus simulating ^^^Bi and ^os^jj decays in the 




(a) Internal wall e^ext ^jqj (MeV) (b) External wall ey-ext Ejoj (MeV) (c) Top and bottom e^ext Ej^j (iUleV) 



Fig. 15. Phase 2 data. The energy sum distribution of e7-external events with in- 
coming 7-rays detected at (a) the internal wall, (b) the external wall and (c) the 
end-caps at the top and bottom of the detector. 




(a) e-crossing Ej^j (MeV) (b) e-j^ext E^q^ (MeV) (c) e7-ext E^ (MeV) 



Fig. 16. Results of the Phase 2 fit for the whole detector: a) energy sum of crossing 
electron events, b) the energy sum of the electron and 7-ray for e7-external events, 
c) the detected electron energy for e7-external events. 



21 



gap between the iron shield and the tracking chamber walls. For radon the 
activity of 11 Bq/m^ is obtained in good agreement with the results of the 
radon monitoring in the LSM. The value of the thoron activity which agrees 
with the NEMO 3 data is 0.22 Bq/m^. The mean activity of ^°Co for Phase 1 
is higher by a factor of 1.3 when compared to Phase 2. This is reasonable 
considering the ^'^Co half-life of Ti/2=5.2 y. The low-energy 7-ray flux for 
Phase 1 was found to be twice as high than in Phase 2. As one can see in 
Fig. 17 the background model fits well the Phase 1 data too. 




(a) e-crossing E^;,^(MeV) (b) ei-ext E^;,^(MeV) (c) e^ext E, (MeV) 

Fig. 17. Results of the fit for the Phase 1 data: a) energy sum of crossing electron 
events, b) the energy sum of the electron and 7-ray for e7-external events, c) the 
detected electron energy for e7-external events. 

The total number of crossing electron and e7-external events and their ob- 
served energy distributions are well reproduced for the whole detector. Nev- 
ertheless, one may expect that the background may vary from one sector to 
another due to possible inhomogeneities of the detector materials. The dis- 
tribution of the number of observed events by sector is not uniform but the 
pattern is repeated by Monte Carlo calculations, see Fig. 18. It indicates that 
the non-uniformity is mainly due to the detector's acceptance. A small num- 
ber of dead PMTs and Geiger cells that vary from sector-to-sector accounts 
for this. The difference between the number of observed and expected crossing 
electron and e7-external events per sector does not exceed 10%. 

3.4 Test of the model with ee^-external events 

This type of event is similar to the e7-external one, where an incoming 7-ray 
deposits some part of its energy in the calorimeter and hits a foil, but here two 
electrons are emitted from the foil due to double Compton scattering in the foil 
or due to a M0ller scattering of a single Compton electron. The probability for 
a 7-ray to produce two electrons in the foil is about three orders of magnitude 
lower than for a single Compton electron. So the statistics for this channel are 
rather poor. Requiring the 7-ray energy deposit > 200 keV the total number 
of events observed for 969 days of data collection is 420 compared to 409 events 



22 




10 12 14 16 18 20 
Sector number 



10 12 14 16 18 20 
Sector number 



Fig. 18. Two distributions of the number of events per sector (left e7-external and 
right crossing electrons). The data are given by black dots and the MC simulations 
are shown with a solid line. 

expected according to the external background model. The distribution of the 
energy sum of the two electrons in these events is shown separately for Phase 1 
and Phase 2 in Fig. 19. The total number of observed 2e~ events produced 
by detected external 7-rays as well as their energy sum distribution are well 
reproduced by MC calculations. 
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Fig. 19. Energy sum distribution of the two electrons in ee7-external events com- 
pared to the results obtained with the external background model. 



4 Radioactivity inside the source foils 



4-1 Measurements of the internal Bi activity 



The internal ^^"^Bi contamination of the source foils is measured with BiPo 
decays through the detection of ea events. The energy loss of alpha particles 
in the foil is significantly larger compared to the energy loss in the gas. As a 
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consequence the a-track length depends on the ^^^Po decay location. There- 
fore the delayed track length distribution was used to measure the internal 
impurities of the source foils. 

As in the case of radon activity measurements it is assumed that the ^^^Bi 
is deposited on the wire surfaces. The possibility of a deposition on the foil 
surface is also considered. Phase 2 data is used to minimize the contribution 
from the radon in the tracking detector. Only events from the source foil (see 
Fig. 5a) are used. The selected events are divided into four groups according 
to the delayed track location with respect to the foil (inner or outer side) and 
to the electron track (on the same side of the foil or on the opposite side). 

An example of delayed track length distributions is shown for the selenium 
foils of the first productioiO], Fig. 20. These are fit with five contributions 
corresponding to the different locations of the ^^^Bi. The five are from inside 
the pf3 material, on the two source surfaces and in the wire surfaces on both 
sides of the foil. 

82c 



a track length (cm), Phase 2, ° Se(l) 




20 20 

(a) a inner foil side, (xe same side (b) a inner foii side, ae diff sides 




data 

inner foil surface 
outer foil surface 
inner side wire surfaces 
outer side wire surfaces 
inside myiar fiim 
inside 33 materiai 



20 20 

(c) a outer foil side, ae same side (d) a outer foil side, oce diff sides 

Fig. 20. Delayed track length distributions for the events of Phase 2 with vertices 
on the selenium source foils is shown by points with error bars. The MC simulations 
are shown with histograms. The results for events with both the a track and the 
fast track (e~) on the same side of the source foil are shown in plots (a) and (c) 
while those having the a track and the electron on different foil sides are shown in 
plots (b) and (d). 




Except for the metallic molybdenum (^°°Mo(m)) and the copper foils, the 
NEMO 3 source foils are made of thin mylar films sandwiching the active 
l3(3 material. The mylar films are from three different productions and their 
activities were measured with HPGe detectors before foil fabrication. The 

^ Selenium foils are from two separate enrichment tasks, which are called ^^Se(I) 
and 82Se(II). 
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Table 7 

Measurements of ^^^Bi activity of the source foils, mBq/kg. Al - results of the 
fit for the foil activity excluding the contribution from mylar. A2 - mylar activity 
measured with HPGe detectors, A3 - total foil activity including mylar calculated 
from Al and A2, A4 - foil activity measured with HPGe. The following notations 
are used for different source foil types: ^'^'^Mo(m) — for the metallic molybdenum, 
^'^*^Mo(c) — for the composite molybdenum, ^^Se(I) — for the selenium sample of 
the first enrichment and ^^Se(II) — for the second. 



Source foil 


Al 


A2 


A3 


A4 


lOOMo(m) 


<0.1 




<0.1 


<0.39 




<0.1; 0.30.0±0.07 


<0.67 


<0.15; 0.27±0.07 


<0.34 


S2Sc{I) 


1.0±0.13 


1.7±0.5 


1.1±0.17 


<4.2 


«2Se(II) 


0.4±0.15 


1.7±0.5 


0.53±0.18 


1.2±0.5 * 


96 Zr 


6.4±2.; 7.8±2. 


<0.67 


5.5±1.7; 6.5±1.7 


<16.7 


IBONd 


2.7±0.4 


3.3±0.5 


2.8±0.4 


<3.3 


i30Te 


<0.1 


3.3±0.5 


0.39±0.06; 0.48±0.06 


<0.67 


Te(nat) 


0.13±0.1 


1.7±0.5 


0.28±0.14 


<0.17 




0.6±0.15; 0.67±0.13 


<1.0 


0.65±0.13; 0.59±0.13 


<1.7 


Cu 


<0.1 




<0.1 


<0.12 



* samples of *'^Se(I) and *^Se(II) having a combined mass of 800 g 
have been measured. The *^Se(II) was not measured separately. 



results of these measurements are given in column A2 of Table 7. After the 
source foil fabrication the foil activity was again measured. Results are shown 
in column A4 of Table 7. For most of the measurements only limits could be 
achieved. 

To make the fit the ^^^Bi activity inside the mylar films covering the foils (if 
any) is fixed to the value measured with the HPGe detector. The measurement 
of the internal ^^^Bi contamination of each NEMO 3 source material, using 
this method, is given in column Al of Table 7. When only a limit on the mylar 
activity is available two results are given. The first is obtained with the limit 
value, the second with the mylar activity set to zero. 

Column A3 shows the total ^^^Bi activity of the source foils calculated from 
the values of columns Al and A2. When only a limit is known for the mylar 
contamination by ^^^Bi there is a systematic uncertainty on the total foil 
activity reflected by the difference between two numbers in the column A3. 
When a limit is obtained on the internal foil material activity Al (the case 
of ^^''Te) the first number in the column A3 corresponds to A1=0, the second 
is for the limit value of Al. There is a good agreement between the results 
obtained from NEMO 3 data and those from the HPGe detectors. 



25 



4-2 "^^^Tl inside the source foils 

The presence of a small quantity of ^oS'p^ from the (Jecay chain inside a 

source foil is the origin of the most troublesome background for the neutrino- 
less double beta decay search. Therefore the radiopurity goals of molybdenum 
and selenium foils were very high [1]. The measurements using HPGe detec- 
tors could not reach the required sensitivity and in most cases only limits 
on 208rp| activity were set. Here the measurements of ^^^Tl foil contamination 
performed with NEMO 3 data are presented. 

Events of 677 ^iid 6777 topologies are used. The event selection criteria are 
similar to those described in section 2.2.1 but in this case the event vertex has 
to be on the source foil. All likely backgrounds have been taken into account. 
The most important one is due to the thoron and radon in the detector volume. 
The two event topologies give consistent results and the activities evaluated 
with Phase 1 and Phase 2 data are in a good agreement. The ^os^l activity in 
the /3(3 source foils and in the copper foils based on the total available statistic 
are shown in Table 8. 

The results of HPGe detector measurements are also shown in Table 8 for 
comparison. With the exception of natural tellurium there exists a good agree- 

Table 8 

Number of observed events (N), signal-to-background ratio (S/B), signal efHciency 
(e) and results of the measurements of ^os-pj activity of the source foils compared 
to the HPGe measurements. 



fSP material 


N 


S/B 


e(%) 


A (mBq/kg) 


^HPGe (mBq/kg) 


lOOMo(m) 


666 


2.4 


1.7 


O.llitO.Ol 


<0.13;<0.1;<0.12* 


iooMo(c) 


1628 


1.7 


1.8 


0.12±0.01 


<0.17 


82 Se (I) 


446 


2.0 


2.0 


0.34±0.05 


<0.670 


82Sc (11) 


507 


3.4 


1.9 


0.44±0.04 


0.4±0.13** 


«Ca 


42 


4.1 


1.4 


1.15±0.22 


<2. 


96 Zr 


158 


7.9 


1.8 


2.77±0.25 


<10.;<5.* 


IBONd 


1002 


39.4 


1.8 


9.32±0.32 


10.±1.7 


130Tc 


448 


1.1 


2.0 


0.23±0.05 


<0.5 




495 


1.9 


1.8 


0.27±0.04 


<0.08 


ii6Cd 


196 


0.7 


1.7 


0.17±0.05 


<0.83;<0.5* 


Cu 


()() 


().() 


1.5 


().0:-!±().01 


<0.033 



* different foil samples have been measured 

** samples of 82Se{I) and ^^gc^jj) having a combined mass of 800 g 
have been measured. The 82Se(n) was not measured separately. 



ment between the two measurements. Due to the large data acquisition time 
NEMO 3 achieves better sensitivity. 
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4-3 Measurement of the internal activity of (3 emitters 



The beta activity of a foil is measured by events with a single electron track. 
The sources of radioactivity considered above are not sufficient to explain the 
single electron data. An additional source of electrons in the foils is required. 
The results of the foil radiopurity measurements performed with HPGe detec- 
tors before the foil installation in NEMO 3 are used to decide upon the list of 
contaminants in a given foil. This list typically includes and isotopes from 
^^'^U and ^ssfj]^ decay chains. The internal ^^^Ra and ^^^Ra activities are fixed 
from the results obtained for ^^^Bi and ^o^xi described above. Their daughters 
are assumed to be in equilibrium. A possible foil surface pollution by ^^°Pb is 
also considered. 



HPGe measurements made for the copper foils have revealed that the copper 
is quite clean. Upper limits of 5 mBq/kg and 8 mBq/kg on the contamination 
by ^^^™Pa and respectively were attained. Here the internal activities of 
these beta emitters are determined by a fit to the single electron energy spectra 
obtained from Phase 2. 




Fig. 21. Single electron energy spectra for copper, Phase 1 and Phase 2 data. 



The results of the fit are shown in Fig. 21 demonstrating a good agreement 
with the data. An activity of 1.5 ± 0.1 mBq/kg is found for ^^^"'Pa and 3.7 ± 
0.1 mBq/kg for '^''K. No pollution on the copper foil surface by ^^''Pb was 
evident. In Fig. 21 the ^^"Bi contribution corresponds to the wire surface 
contamination. 



The same method provides the determination of the contaminant activities in 
the j3l3 source foils. The "single electron" channel is the most appropriate for 
pure beta emitters. For isotopes decaying with copious 7-ray emission such as 
^°Co and ^°^Bi, the e7-internal channel is used. 
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5 Test of the background model 



The highly radiopure copper foils (621 g) occupy one of the 20 sectors in 
NEMO 3. They are used to measure the external background. The internal 67 
and 2e~ channels are used to compare the data with the background model 
described above. 



5. 1 Internal 67 events from the copper foils 



In Fig. 22 some distributions of internal 67 events coming from the copper 
foils are compared with the MC simulation based on the background model. 
One can see that the total number of events and the energy spectra are closely 
reproduced. 



J foils, Phase 1 



Cu toils, Phase 1 



Cu toils. Phase 1 




1 1.5 2 2.5 
et^internal E. 

J toils. Phase 2 



3 3.5 
TOT (MeV) 



1 1.5 2 2.5 3 
ey-internal Eg (MeV) 

Cu toils. Phase 2 



1 1.5 2 2.5 3 
e-^internal E^ (MeV) 

Cu toils. Phase 2 




0.5 1 1.5 2 2.5 3 3.5 
e-y^internal E^qt (MeV) 



0.5 1 1.5 2 2.5 3 
e"(^internal Eg (MeV) 



1 1.5 2 2.5 3 
e-y^internal E (MeV) 



Fig. 22. The distribution of the energy sum of the e~ and 7-ray, the single e^ energy 
and the single 7-ray energy of internal e~7 events from copper foils for Phase 1 and 
Phase 2 data. 



The external background is dominant and represents almost 90% of all internal 
67 events in Phase 2. In Phase 1 the radon produces a noticeable contribution, 
and in the energy sum region below 1 MeV there are slightly more events 
than expected. The difference is less than 2% of the total number of events. In 
Phase 2 the radon contribution is significantly lower and the model describes 
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the data very well even at low energies. A similar problem at low energies for 
the Phase 1 data is also observed in the 2e~ channel and may have the same 
origin. 

The internal background contribution in the 67 channel is negligible so the 
data confirms that the copper foils are clean. 



5.2 Two electron events from the copper foils 

Copper is not a source of double beta decay. That allows one to check the 
validity of the background model in the 2e~ channel where double beta decay 
is searched for. To study this channel the events are selected by requiring two 
reconstructed electron tracks emitted from the foil with the correct curvature 
and a common vertex in the foil. The energy of each electron measured in the 
calorimeter is required to be greater than 200 keV. Each track must hit an 
isolated scintillator block and no additional PMT signals are allowed to avoid 
events with 7-rays. The event is also recognized as internal by the time-of-fiight 
difference of the two electrons. 

In Fig. 23 the distributions of the energy sum of the two electrons, the single 
electron energies and angular correlation of two-electron events coming from 
the copper foils are compared to the prediction of the background model. 
In Phase 2 the total number of events (220 observed with 213 expected for 
788 days of data acquisition) have energy and angular distributions in good 
agreement with the MC predictions. In Phase 1, where the radon activity 
level is six times higher, there are 262 events observed for the 374 days of 
acquisition, which is 15% higher than the MC prediction. The excess of events 
is observed in the energy sum distribution below 1 MeV. The difference is 
relatively small and is apparently due to radon. 

The radon contribution is provided by ^^"^Pb ((5^=1.023 MeV) and ^^^Bi 
((5/3=3.272 MeV) decays simulated according to the model based on [14]. 
The observed discrepancy at low energies seems to indicate some imperfection 
in the model for ^^'^Pb and can be due to the shape of the beta spectra. Shapes 
of beta spectra corresponding to allowed decays were generated, while the beta 
transitions are mainly of the first non unique forbidden type. A deviation from 
the allowed shape for this type of transition may be large [19], [20], however 
it is not known for ^^^Pb and ^^^Bi. 

Apart from the problem at low energies in Phase 1 the background model 
successfully reproduces the 2e~ events produced in the copper foils due to the 
external and internal radioactive sources. 
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Fig. 23. Background model prediction compared to the data for 2e events from 
the copper foils. 

6 Summary 

The methods of background measurements in the double beta decay exper- 
iment NEMO 3 are presented. The background is classified as internal and 
external according to its origin. Internal backgrounds come from the source 
foils. The external background is subdivided into two groups. The first is due 
to radioactive sources inside the tracking volume, and the second results from 
the radioactivity outside of the tracking volume. All these backgrounds were 
estimated from the data with events of various topologies and summarized 
below. 

- The external background component due to the presence of radon in the 
tracking chamber is measured using events with a detected electron accom- 
panied by a delayed a-particle track. This topology results from the beta 
decay of ^^^Bi followed by the alpha decay of ^^^Po. 

- The thoron is measurable with the events of 677 and 6777 topologies. The 
method is based on the detection of events with the signal of a 2.615 MeV 
7-ray, typical of the ^os^j] beta decay. 

- The presence of ^^°Pb on the wires is measured with a single electron 
starting from a wire. 

- The external 7-ray fiux coming from outside of the tracking detector volume 
is measured with a crossing electron and e7-external events. 
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- The internal background due to ^^^Bi and ^os^j decays inside the foil is mea- 
sured in a manner similar to that used for radon and thoron measurements 
but requires the event vertex to be located on the foil. 

- The use of single electron events without detected 7-rays and coming from 
a source foil allows one to measure the activity of beta emitters inside the 
foils. 

Very pure copper foils are used for the study of the external background with 
internal 67 events. The validity of the background model is successfully verified 
with 2e~ events coming from the copper foils. 

It has been demonstrated that with the NEMO 3 detector the backgrounds 
are measurable from the experimental data. In particular the activities of the 
two most troublesome background sources for /3/30z/ decay, ^^'^Bi and ^''^Tl, are 
measured with the adequate precision. 
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